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ABSTRACT 

We report the identification of a source coincident with SN 2009kr in HST 
pre-explosion images. The object appears to be a single point source with an 
intrinsic colour V — 7=1.1 and My = —7.6. If this is a single star it would 
be a yellow supergiant of logL/L ~ 5.1 and a mass of 15+4 M . The spatial 
resolution does not allow us yet to definitively determine if the progenitor object 
is a single star, a binary system, or a compact cluster. We show that the early 
lightcurve is flat, similar to IIP SNe, but that the the spectra are somewhat 
peculiar, displaying unusual P-Cygni profiles. The evolution of the expanding 
ejecta will play an important role in understanding the progenitor object. 

Subject headings: supernovae: general — supernovae: individual (SN 2009kr) — stars: 
evolution — galaxies: individual (NGC 1832) 
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Introduction 



In recent years the hypothesis that red supergiants between 8 and 17 solar masses 
explode as Type II Plateau (IIP) SNe has been confirmed by several direct detections of 
progenitors in pre-explosi on im ages, and by setting upper mass limits from non-detections 
(see, for example, Smartt 



20091 ). Two progenitor candidates which show Luminous Blue 
Variable characterics, either in quiescence or outburst, have been found for SNe with clear 
signa tures of circumstellar interaction, the Ibn 2006jc and Iln 2005gl (Pastorello et al. 



2007; Gal- Yam & Leonard 



20091 ). Despite these successes, the total number of supernovae 



for which we have definite progenitor identifications is still low. SN_2009kr was found in 



the spiral galaxy NGC 1832 on Nov. 6 20 09 by Itagaki (Nakano 120091 ). at an unfiltered 
magnitude of 16.0. Tendulkar et al. (120091 ) obtained a spectrum on Nov. 8, and suggested 
that SN 2009kr showed the features of a Iln supernova. The authors reported a blue, 
featureless continuum, with narrow hydrogen Balmer emission lin es an d a Ha line-width 
of 880 km s _1 . A second spectrum was obtained by Steele et al. (120091 ) on Nov. 9, which 



showed weak P-Cygni absorption in the hydrogen lines. Based on this, and indications that 
the narrow Balmer lines seen by Tendulkar et al. were in fact produced by a nearby HII 
region, Steele et al. claimed SN 2009kr to be a young Type IIP. 

In this paper, we have taken the host galaxy of SN 2009kr to be at a distance of 
26.2 ± 1.8 Mpc (from the recessional velocity, correcting for Virgo-centric infall; with values 
taken from NED), which corresponds to a distance modulus (m — M) of 32.09 ± 0.15 
mag. Estimates from the Tully-Fisher relation howev er, give a significantly higher value 



of (m — M) = 32.61 ± 0.43 mag (Willick et al. Il997f ). This value may be unreliable as 



NGC 1832 h as a fa int absolute magnitude, where the T-F relation has a larger scatter. 



Terry et al. (120021 ) give a low value of (m — M) = 30.76 ± 0.41 mag from the method 



of "sosie" galaxies. These two methods bracket our recessional velocity distance; in the 
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absence of any indication of which method is the most reliable in this case, we used the 
recessional velocity based distance as a compromise, with a conservati ve erro r of ±0.5 mag 



to reflect the uncertainty We have used the standard Schlegel et al. (119981 ) relation for 
Galactic extinction, giving Ay = 0.242 and Ai = 0.142 mag (values taken from NED). Host 
extinction is assumed to be negligible, as our spectra of the supernova appear blue and lack 
prominent Nal D lines. We have used the de-projected galactocentric distance (5.6 kpc) of 
the supernova to estimate the metallicity of the region in which SN 2009kr exploded. We 
first determined a characteristic oxygen abundance for the host g alaxy of 12 + log(0/H) 



8.33 ± 0.24 dex from its absolute B magnitude (Pilyugin e t al. 



oxygen abundance gradient taken from Boissier & Prantzos (120091 ). —0.50 dex / R25, we 



20041 ). Using the radial 



find a metallicity of 12 + log(0/H) = 8.06 ± .24 de x. While this is low compared to the 
location of other Type IIP SNe (Smartt et al. 120091 ). we note that e mpirica l metallicity 



calibrations suffer considerable uncertainties (see for example, Bresolin 



20081 ). 



2. Observations and data analysis 

The host galaxy of SN 2009kr was observed on 2008 Jan 11 (~ 660 days before 
explosion) with the Wide Field and Planetary Camera 2 (WFPC2) onboard the Hubble 
Space Telescope as part of HST Program 10877. Data were reduced and calibrated by 
the on-the-fly calibration pipeline, and downloaded from the MAST archive at STSct. 
Pre-explosion images consisted of two 230s exposures in the F555W filter, and two 350s 
exposures in the F814W filter. Each pair of images was combined with the CRREJ routine 
within iraJ^I to remove cosmic rays. The site of SN 2009kr fell on the WF3 chip, which has 



x http: / / archive.stsci.edu/ 
2 http: / /iraf. noao.edu 
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a pixel scale of 0.1 "/pixel. A 2340s on-source integration in the Ks filter of SN 2009kr was 
obtained with Naos-Conica (NaCo) on the VLT on 2009 Nov. 21 as part of our supernova 
progenitor prograno The S54 camera was used (pixel scale of 0.054") across a 56" x 56" 
field of view and the supernova itself {my ~ 15 mag) was used as a natural guide star to 
provide adaptive optics correction for the image. Data were reduced using standard iraf 
routines. The pre-explosion WFPC2 image and post-explosion NaCo images are shown in 
Figd 

To determine the position of SN 2009kr in the pre-explosion images, 18 sources 
common to both frames were identified. Their centroids were measured with aperture 
photometry and the resulting list of matched coordinates was used to determine a 
geometrical transformation between the pre- and post-explosion images with iraf geomap. 
Translations and independent rotation and scaling in x and y were allowed for. Six outliers 
lying more than one pixel outside the fit were rejected and the fit re-calculated. The final 
rms error in the transformation (as taken from the output of GEOMAP) was 34 mas. The 
position of the supernova was measured using the three different centering algorithms within 
phot (gaussian, centroid and ofilter) with the mean of the three results being taken as the 
position and the standard deviation of the results as the uncertainty (39 mas). The same 
procedure was used to determine the position and uncertainty (3 mas) in the progenitor 
candidate. The position of the supernova as measured in the post-explosion NaCo image 
was then transformed into WFPC2 coordinate system using the geoxytran task and 
the transformation determined previously. The uncertainties in the progenitor candidate 
and supernova positions were added in quadrature together with the uncertainty in the 
transformation to give the total uncertainty in the procedure. The separation between 



3 Based on observations collected at the European Organisation for Astronomical Research 
in the Southern Hemisphere, Chile, Programs 083.D-0131 & 184.D-1140 



- 8- 



the progenitor candidate and the supernova was found to be 6 mas, which is well within 



the total uncertainty of 52 mas. We t 



ius c onclude that the source indicated in Fig. [T] is 



coincident with SN 2009kr. Li et al. (120091 ) identified the same progenitor in the archival 
WFPC2 images, using an alignment to a post-explosion image obtained with the CFHT 
and Mega-Cam. 

To characterise the pro genito r, we carried out PSF-fitting photometry with the 



hstphot package (Dolphin |2000| ). The coincident source was clearly detected in both 
the F814W and F555W images by hstphot at the 10 a level in F555W and 14 a in 
F814W. We ran hstphot twice, first measuring the sky at each pixel from the mean of its 
neighboring pixels, and secondly recalculating the local sky at the location of the progenitor 
using the pixels immediately outside the photometric aperture. Measured magnitudes were 
m v = 24.53 ± 0.11 and mj = 23.47 ± 0.08 for the first run, and m v = 24.71 ± 0.14 and 
mj = 23.48 ± 0.10 for the second, mj is unchanged within the uncertainties, however there 
is a 0.2 magnitude difference in my. We have taken the magnitudes from the second run as 
more reliable, as the recalculated sky value is likely more appropriate for backgrounds that 
vary rapidly over short distances, such as in this case. To give a more reliable uncertainty 
estimate for our photometry, we have added the error in the output of hstphot in 
quadrature with the difference between the magnitudes given by the first and second runs. 
This gives a final progenitor magnitude of my = 24.71 ± 0.23 and mj = 23.48 ± 0.10 from 
hstphot. With the distance modulus (and its uncertainty) and extinction from Sect, [TJ 
this corresponds to an absolute magnitude of My = —7.62 ± 0.55 and Mi = —8.75 ± 0.51, 
and V - / = 1.13 ±0.25. 

As an independent check of the output of hstphot, aperture photometry was 
performed. We applied the same CTE corrections as hstphot; aperture corrections 
were determined from bright isolated sources. The magnitudes obtained from aperture 
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photometry are 0.2 magnitudes brighter than those given by hstphot; we attribute this as 
likely due to extra flux from nearby sources which the PSF-fitting of hstphot can better 
remove. We obtain V — I — 1.03, which agrees with hstphot within the uncertainties. 

The x 2 an d sharpness statistics in the output of hstphot suggest the source is a 
single star-like PSF; we note however that at the distance of NGC 1832 a single WF chip 
pixel corresponds to ~ 1 3 pc. To attempt to constrain the nature of the progenitor further, 



we used ishape (Larsen Il999l ) to characterize the best PSF fit to the progenitor source. 
ishape is intended to determine the true shape and size of sources down to 10% of the 
PSF size, for sources with a S/N of ~ 30 or greater. As our progenitor is detected at a 
much lower significance level than this, we are wary of relying on the results of ishape in 
this case. This caveat notwithstanding, we used the PSF task within the iraf daophot 
package to create a PSF based on eight bright, isolated sources on the WF3 chip in the 
F814W filter. This PSF was then used within ishape to fit the progenitor. The low x 2 of 
the fit indicates that the source is not clearly extended. 



Bastian et al. (120051 ) suggest that sources with My < —8.6 are likely to be clusters. 
While the source detected is ~ 1 magnitude fainter than this, further extinction in the host 
galaxy could increase the true absolute magnitude of the progenitor candidate. In Fig. 
[2] we plot a colour-magnitude diagram (CMD) of 21 sources detected with a S/N of 5 o 
or above in both filter pre-explosion images. We used the fitting statistics in the output 
of hstphot to try to identify extended sources. While hstphot should return x 2 < 1-5 
and —0.3 < sharpness < 0.3 for single, uncrowded sources, values of x 2 < 2.5 and —0.5 < 
sharpness < 0.5 are reasonable. While all sources within 2.5" were within the sharpness 
limits, four sources were found to have 1.5 < x 2 < 2.5, and have been classed as possibly 
extended on this basis, this is also supported by visual inspection of the images. We also 
inspected pre-explosion R-band and Ha images taken with the Wide Field Camera (WFC) 
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on the 2.5m INT on 20 Mar 2005. The continuum subtracted Ha image was aligned to 
the pre-explosion WFPC2 image, and although the pixel scale of WFC (0.33"/pix) is poor 
compared to WFPC2, identifiable sources of Ha emission are visible. Four of the five 
brightest U-band sources are associated with Ha emission, indicating they are young stellar 
clusters. There is no obviously strong Ha emission at the position of SN 2009kr in these 
images. 



3. Discussion 



The V — I colour of our progenitor after correcting for foreground extincti on is 1.13 



±0.25 mag. Taking the intrinsic colours of supergiants from Drilling & Landolt (120001 ). this 

corresponds to a spectral type of G6, with an effective temperature of 4850K, although the 

uncertainty in colour means it could be as early as Gl or as late as K2. For a G6 star, we 

take a V-band bolometric correction of —0.36 ± 0.2 mag from Drilling & Landolt, with the 

uncertainty corresponding to the variation in bolometric corrections across the range of 

possible spectral type. Applying the bolometric correction to the absolute magnitude of the 

progenitor we find the bolometric magnitude to be —8.0 ± 0.6. Using the standard relation 

between bolometric magnitude and luminosity 

LPVi) = Mf,oi - 4.74 
y L Q -2.5 V ; 

we find a progenitor luminosity of log L/L Q = 5.10 ± 0.24. We show this luminosity and 

temperature on a H-R diagram in Fig. [3} together with the stars evolutionary tracks of 



Eldridge & Tout (120041 ) . We have plotted tracks for two different metallicities, Z = 0.004 
and Z = 0.008 (comparable to the SMC and LMC respectively), we note however that the 
final luminosity of each pair of metallicities is extremely close, indicating that the precise 
progenitor metallicity is not a significant source of error in determining the progenitor 
mass. While the progenitor lies closest to the 20M Q track, it is important to remember 
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that a 2OM model will explode when it has an Fe core, and hence its luminosity will be 
higher by about 0.25 dex. It is more appropriate to compare the progenitor luminosity to 
the luminosity reached by the models at their end point, which corresponds to the helium 
core luminosity at the end of core car bon b urning, i.e. the initial mass - final luminosity 



diagram as discussed in Smartt et al. (120091 ) . We hence find a progenitor mass of 15+1 M . 



From Fig. [21 the suggestion that sources brighter than V ~ -8.5 (marked with a 
line) are clusters seems reasonable. It is important to remember, however, that this does 
not mean that all objects fainter than this magnitude are single sources. It is possible 
that the population of blue sources contains several early type supergiants, and that 
the progenitor of SN 2009kr was originally one of these objects, which was transitioning 
between a blue and a red supergiant when it exploded. A further intriguing possibility is 
that, despite the apparent colour, the p rogenitor of this supernova was in fact a luminous 



blue variable (LBV). Smith et al. (120041 ) have suggested that the bi-stability jump observed 
at a temperature of ~ 21000 K, where the stellar wind properties change from a fast 
wind, with a low M to a slow wind with a high M, may also lead to the formation of 
a pseudo-photosphere. If this occurs, then an early B star could appear to be a yellow 
supergiant from its position on the HR diagram. In Smith et al.'s models, the effect is 
stronger for slightly lower masses (10 M ) than that which we find for our progenitor (15 
M Q ). Furthermore, the models used by these authors are for higher luminosities, log L/L Q 
= 5.7, than that of our progenitor. We also fail to see the indications of strong circumstellar 
interaction in the SN spectrum, which we would expect to see for an object which has 
undergone significant episodes of recent mass loss, such as an LBV. 



Smartt et al. (120091 ) have suggested that there may be a missing progenitor population 
of red supergiants (RSG) compared to what one expects from a typical initial mass function 
and Local Group RSG populations. There have been no detections of RSG SN progenitors 
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above log L/L & ~ 5.1, and one possibility is that stars above this traverse back to the 
blue and become Wolf-Rayet stars. As the progenitor of SN 2009kr was of this luminosity, 
is possible that it wa s in t his bluewards moving phase (as was suggested for SN 2008cn 



by Elias-Rosa et al. 



)j3). 



20091 ). However this scenario does not sit comfortably with the 



low metallicity estimate of around 8.1 dex for the SN environment. Some authors have 
suggested that th e class ical LBV stage could even stretch down to luminosities as low as 5.1 



dex (Smith et al. |2009| ). Hence the apparently hotter photosphere than for an RSG could 
be due to a post-RSG LBV stage during which the core collapses much earlier than most 
theoretical models predict. 

The progenitor appears to be a single source, although a cluster cannot be ruled out, 
except by observing the disappeara nce of the progenitor several years hence (e.g. Gal- Yam 



& Leonard 



2009; Maund & Smartt 



20091 ). The pr ogenitor of SN 2004et was originally 



proposed to be a yellow supergiant (Li et al. 



20051 ) of similar luminosity to the progenitor 



detected here. However, Crockett et al. (2009, in preparation) show that this object was 
not a single star, but a blend of three sources. The spatial resolution of the NGC 6946 
images (the galaxy in which SN 2004et exploded, 0.8" at 5.9 Mpc) of 23 pc is similar to 
the resolution of the WFPC2 images. A binary system could also mimic the colour of the 
progenitor candidate, with a massive, early type star as a companion to a RSG which 
exploded. At the distance of NGC 1832, even a wide, non-interacting binary could appear 
as a single source. When the SN fades, future ACS or WF3 images will determine the 
nature of the source. 

A sequence of spectra spanning a period of about 1 month is shown in Fig. [4] (top 
panel). Early-time spectra show a very blue continuum and weak lines of H and He I. Later, 
with the H envelope recombination onset, the continuum becomes redder and metal lines 
(especially Fe II and Sc II) become prominent. In the middle panel a ~30 d spectrum of 
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SN 2009kr is compared with spectra of SNe 1987A, 2005cs and 2004et at a similar phase. 
The spectrum of SN 2009kr shows a relatively red continuum and the classical Fe II lines 
that are ususally visible in Type II SNe during recombination. However, the weakness of 
Ca II lines in this phase, and the two-component absorption profile of Ha are puzzling. The 
latter can be explained with a very prominent high velocity (~ 15000 km s" 1 ) H component 
in addition to the canonical moderate- velocity (~ 7300 km s _1 ) absorption. Alternatively, 
the bluer shoulder can be interpreted as due to line blending (e.g. with N II). Similar 



profiles have bee n alre ady observed in SNe IIP (e.g. in SN 1999em, see Baron et al. 



2000 



Dessart & Hiller 120061 ). but always with the redder component dominating over the bluer. 
The dominance of emission over absorption in Ha is reminiscent of Type II Linear SNe, 
and would support the stripped (to some extent) envelope progenitor. 

The R-band absolute light curve of SN 2009kr spanning a period of ~ 35 days after 
the core-collapse is shown in Fig. H] (Bottom panel), and is compared with the light curves 
of t hree Type II SNe with well-constra ined progen itor information: SN 1987A (Menzies et 



al. 



1983), SN 2005cs (Pastorello et al. 



2006 



20091) and SN 2004et (Maguire et al. 2009). 



The flat light curve of SN 2009kr suggests that this object is similar to Type IIP SNe. Its 
light curve is only marginally brighter than that of SN 2004et and it is about 2 mags more 
luminous than SN 2005cs. The plateau in SNe II P is at tributed to the recombination of the 



19891 ) did not show a proper plateau, but 



cooling H envelope. SN 1987A (e.g. Arnett & Fu 
a delayed and broad light curve maximum attributed to the small progenitor radius at the 
time of the expl osion. The progenitor of SN 1987A was indeed a compact blue supergiant 



(Gilmozzi et al. 119871 ). whilst more canon ical Type IIP SNe are well-known to have red 
supergiant progenitors (e.g. Smartt |2009i ). In this context, a relatively compact yellow 
supergiant progenitor like that of SN 2009kr is expected to produce a plateau in the light 
curve that slowly increases its brightness with time, although without reaching the extreme 
behaviour of the light curve of SN 1987A. This behavior, however, is not observed in the 
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light curve of SN 2009kr until phase ~ 35 days, although further observations are necessary 
to establish if it will continue to evolve like a normal SN IIP. 

This work was conducted as part of a EURYI scheme award (esf.org/euryi). Data were 
obtained from the Multimission Archive at the Space Telescope Science Institute, operated 
by the Association of Universities for Research in Astronomy, Inc., under NASA contract 
NAS5-26555; the NASA/IPAC Extragalactic Database (NED) and the Isaac Newton 
Group Archive maintained by the CASU at the Institute of Astronomy, Cambridge, and 
observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the 
island of La Palma by the Fundacin Galileo Galilei of the INAF (Istituto Nazionale di 
Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de 
Astrofisica de Canarias. We thank J.S. Vink for his suggestion of a pseudo-photosphere as 
a possible explanation for the yellow colour of the progenitor. MF acknowledges funding by 
the Department of Employment and Learning NI, and SM by Academy of Finland (Project 
8120503). KT has received support from Hungarian OTKA Grant K76816 and from the 
Hungarian Eotvos Fellowship. ST acknowledges support by the Transregional Collaborative 
Research Centre TRR 33 "The Dark Universe" of the German Research Foundation (DFG). 

Facilities: HST (WFPC2), VLDYepun (NACO), TNG (Dolores), NOT (ALFOSC), 
CAO:2.2m (CAFOS), NTT (EFOSC) 
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Fig. 1. — Pre- and post- explosion images of the site of SN 2009kr. Each image is centered 
on the supernova location, with the WFPC2 image being transformed to the same scale and 
orientation as the NIRI post-explosion image, as per the text. Cross marks have been added 
to indicate the progenitor location, circle shows location of sources in Fig. [21 Progenitor 
candidate is located at WFPC2 pixel coordinates (432.8, 149.6). 
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Fig. 2. — CMD of all WFPC2 sources detected at the 5 a level or higher in both filters, 
within an arbitrary radius of 2.5" centered on the progenitor location. Sources which have 
indications that they may be extended from either H a or the output of hstphot have been 
marked with a square or cross respectively. The dashed line at My = —8.5 is an indication 
of the threshold beyond which all sources are likely clusters. Progenitor is marked with error 
bars. 
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Fig. 3. — HR diagram showing evolutionary tracks and the progenitor of SN 2009kr. The 
orange tracks correspond to Z=0.008, while the blue tracks have a metallicity of Z=0.004. 
Red dots correspond to observed red supergiants in the LMC. 
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Fig. 4. — Top panel: Time series of spectra of SN 2009kr from Telescopio Nazionale Galileo 
+ Dolores, Nordic Optical Telescope + ALFOSC, Calar Alto 2.2m telescope + CAFOS, 
New Technology Telescope + EFOSC2, showing the evolution of the SNe, and in particular 
the development of the Ha P-Cygni profile. Narrow components are believed to be from a 
nearby HII region. Bottom left panel: Spectrum of the SN at ~30d, compared to other SNe 
at similar epochs, as discussed in text. Bottom right panel: R-band lightcurve of SN 2009kr 
based on photometry obtained with Liverpool Telescope + RatCam, as compared to other 
SNe (r' images from RatCam have been calibrated to Bessell R). SN 2009kr magnitudes are 
relative to discovery epoch (Nov. 6). 



